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ABSTRACT 

We present an analysis of sixteen galaxy clusters, one group and one galaxy drawn from the Chandra X-ray 
Data Archive. These systems possess prominent X-ray surface brightness depressions associated with cavities 
or bubbles that were created by interactions between powerful radio sources and the surrounding hot gas. The 
central galaxies in these systems harbor radio sources with luminosities ranging between ~ 2 x 10 38 -7 x 
10 44 ergs s . The cavities have an average radius of ~ 10 kpc, and they lie at an average projected distance of 
~ 20 kpc from the central galaxy. The minimum energy associated with the cavities ranges from pV ~ 10 55 ergs 
in galaxies, groups, and poor clusters to pV ~ 10 60 ergs in rich clusters. We evaluate the hypothesis that cooling 
in the hot gas can be quenched by energy injected into the surrounding gas by the rising bubbles. We find 
that the instantaneous mechanical luminosities required to offset cooling range between IpV and 2QpV per 
cavity. Nearly half of the systems in this study may have instantaneous mechanical luminosities large enough 
to balance cooling, at least for a short period of time, if the cavities are filled with a relativistic gas. We 
find a trend or upper envelope in the distribution of central X-ray luminosity versus instantaneous mechanical 
luminosity with the sense that the most powerful cavities are found in the most X-ray-luminous systems. Such 
a trend would be expected if many of these systems produce bubbles at a rate that scales in proportion to the 
cooling rate of the surrounding gas. Finally, we use the X-ray cavities to measure the mechanical power of radio 
sources over six decades of radio luminosity, independently of the radio properties themselves. We find that the 
ratio of the instantaneous mechanical (kinetic) luminosity to the 1 .4 GHz synchrotron luminosity ranges from 
a few to roughly a thousand. This wide range implies that the 1 .4 GHz synchrotron luminosity is an unreliable 
gauge of the mechanical power of radio sources. 

Subject headings: galaxies:active - galaxies: clusters: general - X-rays: galaxies - X-rays: galaxies: clusters 



1. INTRODUCTION 

The cooling time of the intracluster gas in the cores of many 
galaxy clusters is shorter than 1 Gyr . In the absence of heat- 
ing, a "cooling flow" (Fabian 1994) is established, in which 
the gas cools below X-ray temperatures and accretes onto 
the central cluster galaxy, where it accumulates in molecular 
clouds and forms stars. Chandra images of cooling flow clus- 
ters have confirmed the existence of inwardly decreasing tem- 
perature gradients and short central cooling times, which are 
the distinguishing characteristics of a cooling flow. However, 
moderate-resolution Chandra and AS CA spectra and high res- 
olution XMM-Newton spectra (e.g. [M akishima et alJ|200il 
iPeterson et"aUl2001tlTamura et aljhoOltlKaastra et alJl2004l) 
do not show the expected signatures of cooling below 2 keV, 
reported to exist in lower resolution data from the Einstein and 
ROSAT observatories. This discrepancy would be difficult to 
understand unless the normal signatures of cooling below 2 
keV are somehow suppressed, or if cooling is indeed occur- 
ring but at rates that are generally factors of 5 - 10 lower than 
expected (e.g. [ M olendi & Pizzolatoll2001t iBohringer et alJ 
I2002L IPeterson et all2003lT Several scenarios have been sug- 
gested that may suppress the cooling flux at low energies, and 
yet maintain the large cooling rates. The list includes differ- 
ential absorption, efficient mixing, or inhomogeneous metal- 



licity distributions (e.g. Fabi an et"aT]l200 1). However, these 
scenarios, as clever as they are, may require fine-tuning and 
are otherwise difficult to prove observationally. In any case, 
they leave the question of the repository for most of the cool- 
ing gas unanswered. 

The more appealing interpretation, which we address in 
this paper, posits that radiation losses are being balanced, 
or nearly so, by heating, implying that the large cool- 
ing rates of the last decade were indeed overestimated. 
This suggestion has its own difficulties. Maintaining gas 
with a cooling time approaching 100 Myr at keV tem- 
peratures almost certainly requires the existence of one 
or more heating mechanisms operating in a self-regulating 
feedback loop. One such mechanism, thermal conduc- 
tion from the hot outer layers of cluste rs, may be ener- 
getically feasible, in some instances ( e.g. | Tuc ker ffi: Rosner 
1 19831: iBertschinger & Meiksinl 119861 IZ akamska & Naravanl 
l2003[ iVoigtetalJ 120021 iVoigt & Fabianl 120041 How- 
ever, it generally requires fine tun ing and can be unstable 
(Bregman & Davidll 988: Soke J2003T) . Moreover, conduction 
operati ng alone at even the Spitzer rate ca nnot offset radiation 
losses ( Voigt et al. 2002; Wise et al. 2004) in all clusters, and 
is therefore unlikely to provide a general solution to the heat- 
ing problem (recent simulations by Dolag et al. 2004 support 
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TABLE 1 
Radio properties. 







a (ref) a 


Sl400 




-Pi 400 


^-radio 


System 


z 


(km s" 1 ) 


(Jy) 


a (ref) b 


(10 24 \VHz~ 1 ) 


(10 42 ergs s" 1 ) 


Cygnus A 


0.056 




1598±41 


0.7 (4) 


11800±300 


700±20 


Hydra A 


0.052 


322±20 (8) 


40.8±1.3 


0.92 (4,13) 


261±8 


19.9±0.6 


A2597 


0.085 


224±19(16) 


1.875±0.056 


1.35 (11) 


35±1 


6.8±0.2 


MKW3S 


0.045 




115.0±3.9 


2.3 (4,13,17) 


0.58±0.02 


3.6±0.1 


A2052 


0.035 


253±12 (19) 


5.50±0.21 


1.2 (1,4,13) 


15.7±0.6 


2.08±0.08 


A133 


0.060 




0.167±0.006 


1.9 (14,15) 


1.54±0.06 


1.94±0.07 


A4059 


0.048 


296±49 (6) 


1.284±0.043 


1.43 (13,18,21) 


6.96±0.2 


1.73±0.06 


A2199 


0.030 


295±6 (7) 


3.58±0.12 


1.37(1,4,17) 


7.43±0.2 


1.56±0.05 


Perseus 


0.018 


246±10(8) 


22.83±0.68 


1.0 (10) 


16.5±0.5 


1.42±0.04 


RBS 797 


0.350 




0.0217±0.0008 




9.0±0.3 


0.78±0.03 


A 1795 


0.063 


297±12 (3) 


0.925±0.028 


0.98(1,4,13) 


9.0±0.3 


0.75±0.02 


M87 


0.0042 


330±5 (2) 


138.5±4.9 


0.81 (4,9,13) 


5.48±0.2 


0.36±O.Ol 


Centaurus 


0.011 


256±11 (5) 


3.8 


0.7 (4,9,14,20) 


1.1 


0.060 


A478 


0.081 




0.0369±0.0015 




0.60±0.02 


0.052±0.002 


M84 


0.0035 


278±4 (8) 


6.00±0.15 


0.63 (4,9,17) 


0.162±0.004 


0.0089±0.0002 


2A 0335+096 


0.035 




0.0367±0.0018 


0.9 (12) 


0.104±0.005 


0.0077±0.0004 


A262 


0.016 




0.0657±0.0023 


0.6 (4) 


0.039±0.001 


0.002 10±0.()0007 


HCG 62 


0.014 




0.0049±0.0005 




0.0021±0.0002 


0.00018±0.00002 



REFERENCES. — (1) Becker, White, & Edwards 1991; (2) Bender, Saglia, & Gerhard 1994; (3) Blakeslee & Tonry 1992; (4) Burbidge & 
Crowne 1979; (5) Carollo, Danzinger, & Buson 1993; (6) Carter et al. 1985; (7) Fisher, Illingworth, & Franx 1995; (8) Heckman et al. 1985; (9) 
Kuehr et al. 1981; (10) Pedlar et al. 1990; (11) Sarazin et al. 1995; (12) Sarazin, Baum, & O'Dea 1995; (13) Slee 1995; (14) Slee & Siegman 
1988; (15) Slee et al. 2001; (16) Smith, Heckman, & Illingworth 1990; (17) Spindrad et al. 1985; (18) Taylor, Barton, & Ge 1994; (19) Tonry 
1985; (20) Wright et al. 1994; (21) Wright et al. 1996 

a When no velocity dispersion was available, (a) = 280 km s" 1 was adopted. References are in parentheses. 

^The spectral index is defined so that S ~ v~ a . When no spectral index was available, a = 1 was adopted. References are in parentheses. 



this conclusion). Add itional he at sources, such as c osmic rays 
( Bohri nger & Morfill 1988llLoewenstein et alJl99lT) . and su- 
pernova explosions ( McNamar a"etai*l2004l) may contribute to 
heating the gas. Nevertheless, these mechanisms are generally 
incapable of balancing radiation losses. 

In this paper, we evaluate whether the mechanical 
energy generated by active galactic nuclei (AGN) can 
balance radiation losses in cluster cores. This pos- 
sibility, wh ich has a substantial legacy in the litera- 
ture (e.g. iTab or & Binnev 1993; Binnev & Tabor J19951 
Tucker & David) 119971 iCiotti & Ostrikerl 1200 It iSoker et all 
1200 II) . has been rejuvenated by the crisp, new Chandra im- 
ages of clusters showing the keV gas being displaced by ra- 
dio sources harbored by central cluster galaxies. The now 
ubiquitous signature of these interactions are X-ray surface 
brightness depressions projected on the r adio lobe emis- 
sion at 1.4 GHz, as is seen in Perseus jBohringer et alJ 
[T993t IFabian et al J l2000l iSchmidt et all 120021 IFabian et all 
l2002bl l2003alH). Cvgnus A (iCarilli et all 119941). and Hy- 
dra A ( IMcNamara et all2000tlDavid et al.l2001tlNulsen et al.l 
120021) . The displacement of the gas creates a low-density, ris- 
ing bubble in pressure balance with the surrounding medium. 
X-ray surface brightness depressions that have no obvious as- 
sociation with the bright radio emission at 1 .4 GHz, the so- 
called ghos t cavities, have also bee n found, such as those in 
Abell 2597 ( McNamara et al. 2001 ) and the outer depressions 
in Perseus (Fabi an et alJ20 00). These depressions are thought 
to have been created by interactions that occurred in the more 
distant past, but whose radio emission has faded over time. 

This general scenario has been modeled theoretically 
using a variety of hydrodynamical, magnetohydrody- 
namical, and analytical techniques, which have suc- 
cessfully reprod uced the gross charac te ristics of the 
cavities (e.g. iGull & Northoverl 119731 IChurazov et al.l 



pOOlt BrUggen & Kaiserl l200lt lOuilisetalJ 2001; 
Brighenti & Matthews! 120021; IBriiggen et alJ 2002 : 
Reynolds, Heinz, & Begelman 2002; Brugg enl 12003 ; 
De Yound 120031: iBasson & Alexander! 120031 [Binnevl 12004 
Kaiser & Binnevl 120031 iMathews et alJ 120031; lOmma et all 
2004; Robins on et alJl2004l) . The models do not, however, 
predict with any certainty the amount of mechanical energy 
provided by radio sources of a given luminosity, nor their 
frequency of recurrence. Wherea s the cavities in sev eral 
individual objec ts, such as P erseus ( Fabian et al. 2003a|) and 
Hydra A (McNamar a et alJ 12000), contain enough enthalpy 
to balance cooling, at least for a short period of time, a 
systematic survey of cavities in systems with a broad range of 
properties is required to determine whether this is generally 
true. In this paper we address this question by setting 
observational limits on the energetics and ages of cavities 
in 18 systems taken from the Chandra Data Archive. We 
adopt H = 70 km s" 1 Mpc" 1 , Q M = 0.3, and £1 A = 0.7 in all 
calculations throughout this paper. 

2. THE SAMPLE 

Approximately 80 systems from the Chandra Data Archive 
were visually inspected for surface brightness depressions. Of 
these, we selected the 18 systems having well defined surface 
brightness depressions associated with their radio sources (see 
Tabled. Of the 18 systems, 16 were imaged with the ACIS- 
S3 detector and two with the ACIS-I3 detector (RBS 797 
and MKW 3S), with exposure times ranging from 12 ksec 
(RBS 797) to 50 ksec (HCG 62). The sample consists of 16 
galaxy clusters, one galaxy group (HCG 62), and one giant 
elliptical galaxy (M84), ranging in redshift from z = 0.0035 
(M84) to z = 0.35 (RBS 797), and in X-ray luminosity from 
- 10 41 ergs s" 1 (M84) to 4 x 10 45 ergs s" 1 (RBS 797). We 
avoided depressions with questionable association with a ra- 
dio source, and we excluded clusters in which there is clear 
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evidence for merging, since merging clusters often show com- 
plex structure that can be mistaken for a radio-induced cavity. 
All clusters in our sample were previously reported in the lit- 
erature as containing cavities likely to be associated with radio 
bubbles. 

3. DATA REDUCTION AND ANALYSIS 

3.1. Radio Analysis 

The systems in our sample have a wide range of radio prop- 
erties, from powerful double-lobed FR II radio sources with 
luminosities of ~ 7 x 10 44 ergs sec" 1 (Cygnus A) to weak 
sources with luminosities of ~ 2 x 10 38 ergs sec -1 (HCG 62). 
Table [2 gives the radio properties of our sample. The ra- 
dio power at v = 1400 MHz was calculated as P v = 4irD^Su, 
where D v = D L (l+z)~ (l+a)/2 ivon Hoernerl 19741) . The total ra- 
dio luminosity was calculated by integrating the flux between 
vi = 10 MHz and v 2 = 5000 MHz as 

L rad = 47rD 2 S„ / (y/v Q ) a dv, (1) 

J U\ 

where we have assumed a power law spectrum (JS V ~ v~ a , 
where a is the spectral index). We used for S„ the 1400 MHz 
flux from the NRA O VLA 1 Sky Survey (NVSS) catalog 
llCondon et alJ 199 8k except in the case of the Centaurus clus- 
ter where no NVSS data were available. In this case, we used 
the 141 MHz flux from th e Parkes Radio Sources Catalogue 
(Wright & Otrupcek 1990). Spectral indices were taken from 
the catalogs referenced in Table ^ As the derived spectral 
index can vary depending on the frequencies used, we have 
adopted a weighted average of the available spectral indices. 
In cases in which the spectral index was not available, a value 
of a = 1 was adopted. 

3.2. X-ray Analysis 

The X-ray data were obtained through the Chandra Data 
Archive and were reprocessed with CIAO, version 2.3, using 
CALDB, version 2.21. The charge transfer inefficiency (CTI) 
correction was applied during reprocessing of the level 1 event 
file. Blank-sky background files were used for background 
subtraction for all clusters. 2 The background files were nor- 
malized to the count rate of the source image in the 10-12 keV 
band, after the removal of all bright emission. The required 
adjustment was less than 12% for all clusters except Centau- 
rus and Perseus, which both required background adjustments 
of ~ 30%. Spectra with at least 2000 counts were extracted 
in circular annuli centered on the X-ray centroid of the clus- 
ter. Response files were made using the CIAO tools mkrmf 
and mkwarf. We attempted to correct the resulting ARFs for 
the quantum efficiency degradation problem using the corrarf 
tool. 3 However, upon spectral fitting, we found that ~ 75% 
of our sample was overcorrected by corrarf. Therefore, we 
present our results without th e correction applied ( for a dis- 
cussion of this problem, see Voigt & Fabian 2004). In gen- 
eral, the largest effect of using corrarf was on the cooling 
rates, which increased on average by a factor of 2 after the 
correction was applied. 

1 The VLA (Very Large Arry) is a facility of the National Radio Astron- 
omy Observatory (NRAO). The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universi- 
ties, Inc. 

2 See http://asc.harvard.edu/contrib/maxim/acisbg/ 

3 See http://asc.harvard.edu/cal/Acis/Cal_prods/qeDeg/ 



To find radial temperatures and densities, we deprojected 
the spectra extracted above. The deprojection was performed 
assuming spherical symmetry and using the PROJCT model 
in XSPEC 11.2.0 with a single-temperature plasma model 
(MEKAL) and foreground absorption (WABS), fitted be- 
tween energies of 0.5 and 7.0 keV. The foreground column 
density, A^h , was tied between annuli and allowed to vary. The 
MEKAL abundance was also free to vary. The redshift was 
fixed to the value given in Tabled The density was then cal- 
culated from the normalization of the MEKAL component, 
assuming n e = 1 .2«h (for a fully ionized gas with hydrogen 
and helium mass fractions of X = 0.7 and Y = 0.28). The pres- 
sure in each annulus was calculated as P = nkT, where we 
have assumed an ideal gas and n rj 2n e . 

The luminosity of the cooling gas inside the cooling ra- 
dius is needed to investigate whether or not AGN heating 
can balance cooling. Within the cooling radius, radiative 
energy losses must be replaced to prevent the deposition of 
large quantities of cool gas. We define the cooling radius as 
the radius within which the gas has a cooling time less than 
7.7 x 10 9 yr, the look-back time to z = 1 for our adopted cos- 
mology. This redshift is roughly the distance to which clus- 
ters have been found with properties similar to present day 
clusters. The corresponding lookback time should then ap- 
proximate the time a cooling flow has had to establish itself. 
The cooling time was cal culated using the cooling curves of 
Bohringer & Hensler ( 1989). Table|2]gives the values of r coo i 
for each cluster. Within this radius, we performed deprojec- 
tions by fitting both a cooling model and a single temperature 
model to the spectra. 

In order to obtain a spectroscopic estimate of the cooling 
luminosity, we performed the deprojection using a cooling 
flow model (PROJCT*WABS*[MEKAL+MKCFLOW]), fit 
between 0.5 and 7.0 keV. To force all cooling to be within 
the cooling radius, the MKCFLOW model was used only 
inside the cooling radius and was set to zero outside. The 
MKCFLOW low temperature was fixed to 0. 1 keV, resulting 
in an estimate of the luminosity of gas cooling to low temper- 
atures. Within each annulus, the MEKAL and MKCFLOW 
abundances were tied together, and the MKCFLOW high tem- 
perature was tied to the temperature of the MEKAL com- 
ponent. Lastly, the column density (A^h) was tied between 
annuli and allowed to vary. The spectroscopic estimate of 
the bolometric cooling luminosity inside the cooling radius, 
L spec , was then calculated from the unabsorbed fluxes ob- 
tained from the MKCFLOW model, integrated between en- 
ergies of ~ 0. 1 - 100 kev. 

To find the total luminosity inside the cooling radius, we 
performed the deprojection using a single-temperature model 
(PROJCT* WABS *MEKAL) fit between 0.5 and 7.0 keV to 
the same spectra used with the cooling flow model, again with 
A^h tied between regions and allowed to vary. The unabsorbed 
fluxes from the MEKAL components for the annuli within the 
cooling radius, extrapolated between ~ 0.1— 100 kev, were 
used to find the bolometric luminosity of the X-ray emitting 
gas, Lx- Table|2]gives Lx and L spec (with la errors estimated 
by XSPEC) for each object in our sample. Typically, L spec is 
approximately 10% of L x f° r our sample. Our values for L spec 
and Lx are in reasonable agreement with published values for 
most of our sample. 
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TABLE 2 

Results of X-ray spectral modeling. 







PROJCT*WABS*MEKAL 


PROJCT*WABS*(MEKAL+MKCFLOW) 




'cool 


£*(< ''cool) 


M L spc c« ^cool) 


System 


(kpc) 


X 2 /dof (10 42 ergss-') 


X 2 /dof M yr 1 (10 42 ergss-') 



RBS 797 
A478 


191 
150 


154/183 
3704/2796 


4500+™ 
1220±60 


151/180 
3671/2792 


880™ 
150+!$ 


1 200+ 4600 
1 -1100 

180+g 


Perseus 


102 


9348/5717 


670 ^o 


9165/5703 






A1795 


137 


2476/1721 


490±30 


2487/1718 


18+12 


11+9 


A2597 


129 


1341/1149 


430+4JJ 
410±30 


1329/1145 


59 ±40 


28±19 


Cygnus A 


78 


3062/2234 


3062/2230 


< 8 


<6 


2A 0335+096 


122 


2439/1903 


290±20 


2357/1897 


120 ±30 


44+12 

-u 


Hydra A 


100 


1658/1486 


250±15 


1652/1482 


14+?, 


A2199 


113 


2422/1971 


150±10 


2421/1966 


2 + ™ 

—1 .9 

25 ±6 


1 2 +4: 5 


A133 


92 


1624/1160 


95±5 


1598/1158 


11±3 


MKW3S 


88 


2242/2037 


92±7 


2242/2033 


<2 


< 1 


A2052 


101 


2934/2035 


84+« 


2867/2030 


12+4 


4 4+1-7 


A4059 


99 


1396/1109 


7lJ 


1296/988 


6.r«? 


4 7 +3 - 9 


Centaurus 
A262 
M87 
HCG 62 


62 
66 
35 
33 


5100/2177 
1978/1169 
15540/5714 
1008/626 


30±2 
13.8±1.0 

9 8 +0 - 8 
2.0±0.2 


4932/2171 
1923/1165 
15622/5688 
961/624 


10.2*" 

1 8+1-2 

\> 
1 - 1 -0.2 


1 - J -0.4 

62 m34 
-Q18 

20 + * 04 

-Q-Q3 


M84 


14 


682/450 


0.09±0.01 


701/448 


0.06 ±0.05 


uul '-0.004 



4. X-RAY SURFACE BRIGHTNESS DEPRESSIONS 

In total, 36 surface brightness depressions or cavities were 
identified in the 18 systems. Table[3]lists the cavity properties. 
The cavities are classified as either radio-filled or radio-faint 
ghosts, depending on the presence of 1400 MHz or higher 
frequency radio emission inside the cavity. The cavities are 
classified as radio-filled if there is a direct anticorrelation be- 
tween the radio emission at 1400 MHz or higher and the X- 
ray emission, such that the radio emission fills preferentially 
the X-ray surface brightness depressions. The cavities clas- 
sified as ghost cavities, while possibly possessing significant 
radio emission at frequencies at or below 1400 MHz, do not 
show the anti-correlation between the high-frequency radio 
emission and the X-ray emission. Our classification scheme 
relies heavily on the availability of high-resolution radio im- 
ages at several frequencies. However, the radio data available 
are inhomogeneous, and classifying the objects lacking high 
-resolution radio images was challenging. The poor radio im- 
ages available for A262, RBS 797, and HCG 62, in particular, 
led us to classify them as ghosts (see E. L. Blanton et al. 2004, 
in preperation). 

For each cavity, a size and position were measured, assum- 
ing the cavity extends to the inner edge of any bright sur- 
rounding emission. The projected shapes of the cavities were 
measured by eye as circles or ellipses from the exposure- 
corrected, unsmoothed images. This is a qualitative mea- 
surement, the accuracy of which depends on the signal-to- 
noise ratio of the image and on the contrast of the cavity with 
its surroundings. To distinguish between the poorly defined 
and well-defined cavities, we have assigned a figure of merit 
(FOM) to each cavity ranging from 1 for the best-defined cav- 
ities - those with surrounding bright rims - to 3 for the worst- 
defined ones without bright rims. 

In the analysis that follows, we assumed that the cavi- 
ties are bubbles devoid of gas at the local ambie nt tempera- 
ture JMcNamara et alJl2000t iBlanton et alll2003l) . Their vol- 
umes were calculated assuming spherical or prolate ellip- 



soidal shapes, with semimajor axis a and semiminor axis b. 
The errors in the volumes due to projection were estimated by 
allowing each bubble to have an intrinsic a jb as large as that 
of the most eccentric cavity observed in the sample, {a / Z?) 
The upper and lower limits are calculated assuming either 
oblate or prolate symmetry. In this sense, spherical bubbles 
have the greatest range of possible volumes, while projected 

7 /7 

ellipses with an a/£> = (a/^) max have the smallest range. The 
pressure and temperature of the gas surrounding the cavity 
were taken to be the azimuthally averaged values at the pro- 
jected radius of its center. The work done on the surrounding 
medium by the cavity is then simply Wt, u b = pV, if it expands 
slowly compared to the sound speed. 



4.1. Cavity Ages 

The age of each cavity was calculated in three ways. First, 
as the time required for the cavity to rise the projected dis- 
tance from the radio core to its present location at the speed 
of sound, v Cs = (-fkT / iimn) [ ^ 2 , where we have taken 7 = 5/3 
and [i = 0.62. The cavity age is then 

t Cs = R/v c , = R^/nmn/^kT, (2) 

where R is the projected distance from the center of the bub- 
ble to the radio core. Thi s scenario is f avored in the computa- 
tional modeling of Omma et al] J2004I) . in which the bubble is 
produced by a high-momentum jet from the AGN instead of 
rising buoyantly. Second, the age was calculated as the time 
required for the cavity to rise buoyantly (bubble-like) at its ter- 
minal velocity v t ~ (2gV /5C) 1//2 , where V is the volume of the 
bubble, S is the c ross section of the b ubble, and C = 0.75 is the 
drag coefficient dChurazov et a l. 2001). The gravitational ac- 
celeration was calculated using the stellar velocity dispersion 
of the central galaxy, under the approxi mation that the g alaxy 
is an isothermal sphere, as g ps 2a 2 /R (Binnev & Tremaine 
119871) . Published values of the velocity dispersion were used 
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TABLE 3 
Cavity properties. 





Cavity 


a h 


If 


R i 


pv 


i 1 7 vrl 
yiv yi) 


^buoy 






Oy aLClll 




(Vnr\ 
IAPW 


l K PV 


(\cr\r\ 


(10 57 ergs) 




Tin 7 vri 


ref 


RBS 797 


G(2) 


9.7 


9.7 


19.5 


190j|f° 
190!™ 


1.4 


3.6 


6.8 


18 




G(2) 


13.4 


8.5 


23.8 


1.7 


5.2 


7.8 




A478 
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a Radio filled cavities are denoted by "F," radio-faint ghosts are denoted by "G." The FOM gives a relative measure of the cavity's contrast to 
its surroundings: (1) high contrast: bright rim surrounds cavity; (2) medium contrast: bright rim partially surrounds cavity; and (3) low contrast: 
no rim or faint rim surrounds cavity. 

^Projected semimajor axis of the cavity. 

c Projected semiminor axis of the cavity. 

^Projected distance from the cavity center to the radio core. 

e The errors in pV include an estimate of the projection effects; see the text for details. 



when available (see Table otherwise, the average value 
((a) = 280 km s" 1 ) was adopted. The cavity age is then given 

by 

t huoy = R/v t ~R^SC/2gV. (3) 

Finally, the age was calculated as the time required to 
refill the displaced volume as the bub ble rises upward 
aMcNamara et alJ2000tlNulsen et all2002l) : 

t r ~ 2Ry/r/GM(R) = (4) 

where r is the radius of the cavity [for ellipsoidal cavities, 
r=(ab) x l\ 



In general, the ages calculated using the speed of sound are 
the shortest, those based on the refilling time scale are the 
greatest, and those calculated using the terminal velocity lie in 
between. The instantaneous mechanical luminosity per cavity 
or cavity pair is then L mec h = Whub/t, where t is the age of the 
bubble. 

This approximation to the mechanical luminosity of the ris- 
ing bubbles is highly uncertain. In calculating L mec h> a mea- 
surement of the average timescale for the radio source to re- 
plenish itself is required. This timescale is unknown for each 
source, but perhaps this can best be measured in objects with 
clearly defined ghost cavities and a detached central source, 
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FIG. 1. — Left: Mechanical luminosity vs. radio luminosity. The symbols and wide error bars denote the values of the mechanical luminosity calculated using 
the buoyancy timescale. The short and medium-width error bars denote upper and lower limits of the mechanical luminosity calculated using the sound speed and 
refill timescales, respectively. The different symbols indicate FOMs of {circles), 2 (triangles), and 3 (squares). Filled symbols denote radio-filled cavities, and 
empty symbols denote ghost cavities. Each point represents the sum of mechanical luminosities of each bubble type. The best-fit lines are shown for the entire 
sample (dashed line) and for the radio-filled cavities only (dotted line). Right: The ratio of mechanical luminosity to uP v at 1400 MHz vs. radio luminosity for 
the radio-filled cavities only. 



where the duty cycle is clearly evi dent. This situation has 
been n oted in two objects, Pe rseus ( Fabian et alJl200(J) and 
A2597 ( McNa mara et aU2 00 1 ) whose ghost cavities range in 
age between ~ 5-8 x 10 7 yr, with a likely age of ~ 10 8 yr 
when projection is taken into account. Furthermore, of the 80 
or so clusters we searched for cavities, 16 were found to have 
them. If one assumes, for the moment, that all central cluster 
galaxies produce bubbles at a similar rate, then we are seeing 
clusters in an on state only ~ 20% of the time. Therefore, the 
average time elapsed between the production of bubble pairs 
could be as large as 5 x 10 s yr, at least in some clusters. Our 
search did not discriminate between cooling and noncooling 
clusters, and it is likely, although it has not been proved, that 
bubbles are produced preferentially in cooling flows. If true, 
our estimate of ~ 10 8 yr between outbursts is probably closer 
to the truth for most cooling flow clusters. Then the time aver- 
aged mechanical luminosity discussed below has been overes- 
timated in most objects by factors of 2-3. On the othe r hand, 
the 333 MHz radio map of Hydra A JLane et all20 04) shows 
2' plumes extending several times further from the AGN than 
the X-ray cavities (and the 4 GHz radio image). It also shows 
an outer radio lobe 4' north of the AGN that coincides with 
a feature in the X-ray image that hints of a distant cavity. In- 
terpreted as above, the outer feature would give a very long 
interval (~ 10 9 yr) between bubbles, but the plume indicates 
that such an outer bubble has been followed (perhaps some 
time later) by an extended period of continuous radio activ- 
ity. In that case the cavities may just be the latest in a series, 
but we are failing to detect most of the remnants. This would 
make our estimate of L mec h closer to the truth. 

5. DISCUSSION 

5.1. Trends with Radio Luminosity 

In Figure n we present two plots showing the mechanical 
luminosity versus the total radio luminosity (left) and the total 
radio luminosity versus the ratio of the mechanical luminosity 



to the monochromatic, 1.4 GHz radio luminosity (right). In 
each plot we distinguish between radio-filled and ghost cav- 
ities, shown with filled symbols and empty symbols respec- 
tively. The "error bars" for each point reflect the range of 
instantaneous mechanical luminosity implied by the range in 
possible ages. The data are taken or are derived from Tables^ 
and|3] 

The left-hand panel of Figure ^ shows a trend between the 
radio luminosity and mechanical luminosity, with the sense 
that more luminous radio sources tend toward larger mechan- 
ical luminosities. This trend seems to be shared by both the 
radio-filled cavities and the ghost cavities, in spite of the use 
of the current central radio power for both the filled cavities 
and the ghosts, to which the current central source may be 
unrelated. No segregation by FOM is seen. The relation be- 
tween the two luminosities appears to be roughly a power law. 
To quantify this relation, we used a linear least-squares fit to 
the logarithms of the data, with errors in mechanical luminos- 
ity given by the extreme values for each system. We show in 
Figure H the best-fit line for the entire sample (dashed line), 
given by 

t -^(^ 25±3 (r \0.44±o.06 ,~ 



and for the radio-filled cavities only (dotted line), given by 

, 1f vl8±4/r l0.6±0.1 ,f-. 

Anech - 1U (/-radio j ■ \S>) 

In both cases, the mechanical luminosity scales as the ra- 
dio luminosity to approximately the one-half power over six 
decades of radio power, albeit with large scatter. 

The relative contribution of cosmic scatter and observa- 
tional uncertainty is hard to judge without precision radio 
data at a variety of wavelengths and without a better under- 
standing of the bubble production timescale. Nevertheless, 
the existence of this trend demonstrates quantitatively that 
the radio sources are indeed creating the cavities. The radio 
sources are not simply filling preexisting voids in the intra- 
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cluster medium (ICM) created by other processes. Further- 
more, the synchrotron luminosity and mechanical luminosity 
do not scale in direct proportion to each other. This relation- 
ship implies that the synchrotron luminosity cannot be used to 
infer the mechanical power of a radio jet in a simple fashion. 

An important and poorly understood aspect of radio source 
physics is the degree of coupling between the mechanical (ki- 
netic) luminosity of radio sources and their synchrotron lu- 
minosity. This coupling is theoretically tied to the magnetic 
field strength and a ge of the source (see iDe Young||1993t 
iBicknell et al.ll 19971) . neither of which can be measured re- 
liably from radio data alone. Radio sources are inefficient ra- 
diators. The ratio of mechanical power to radio power is typ- 
ically assumed to range between 10 and 100, almost e ntirely 
on theoretical considerations (De Young 1993; Bic knell et al. 
119971) . On the other hand, measurements of the the X-ray 
cavity sizes and surrounding gas pressures provide unique es- 
timates of their ages and mechanical luminosities, indepen- 
dently of the radio properties themselves. We evaluate the 
ratio of mechanical energy to radio power by plotting the ra- 
tio of mechanical power in the bubbles to monochromatic, 
1.4 GHz synchrotron luminosity, assuming IpV of energy per 
radio lobe, against the radio luminosity in the right-hand panel 
of Figure^ This ratio ranges from a few to a few hundred for 
the powerful sources, which is broadly consistent with theo - 
retical estimates (see De Young 1993; Bicknell et al. 1997). 
On the other hand, Abell 478 has a ratio exceeding a few thou- 
sand. To the extent that X-ray cavities provide a good measure 
of the mechanical energy of radio sources, the large variation 
in this ratio indicates that radio luminosity is not necessarily 
a reliable probe of the available mechanical energy. 

There are several factors that can introduce scatter into our 
estimate of the ratio of radio to kinetic power. The most 
important is probably intrinsic differences between the radio 
sources themselves, a consequence of dramatic changes in 
radio luminosity with time. Certainly, if radio outbursts are 
to compensate for radiative losses in cooling flows, then the 
absence of radio emission from some systems requires large 
variations of radio luminosity with time. On the other hand, 
the pV energy of the bubbles alone would tend to underesti- 
mate the the mechanical luminosity of radio sources by fac- 
tors of several if energy dissipating shocks are generated, or 
if the bubbles expand non-adiabiatically (they leak), or if the 
internal energy of the bubbles is boosted with a relativistic 
plasma. 

5.2. Heating by Radio-Induced Cavities 

IChurazov et a noted the conversion of enthalpy of 

the rising bubble into other forms of energy in the cluster at- 
mosphere. Here it is shown that, for an adiabatic bubble, this 
energy is dissipated in its wake. If the mass in the bubble is 
negligible compared to the mass of the gas it displaces, then a 
bubble rises, because the gas falls in around it to fill the space 
it occupied. This process is driven by the potential energy 
released as the surrounding gas moves inward. The energy 
is first converted to gas kinetic energy, then diss ipate d in the 
wake of the rising bubble. In the notation of § 14.11 the po- 
tential energy released when the bubble rises a small distance, 
dR, is 



is the gas pressure. This gives a differential equation for the 
energy dissipated in the bubble wake 



SW = pVgSR = -V^r-SR, 
dR 



(7) 



dW 
~dR 



- V ± 
dR' 



(8) 



If the bubble is adiabatic, with ratio of specific heats 7, then 
pyi = constant and this equation can be integrated to give 
the energy dissipated as the bubble rises over a large distance, 
from Rq to R\, 



AW = —!-(j>oV -piVi)=H -Hi. 



7" 



(9) 



where p is the gas density, and we have used the equation 
of hydrostatic equilibrium to replace pg = -dp/dR, where p 



Here, subscripts and 1 label quantities at the corresponding 
radii, and the enthalpy of the bubble is H = jpV /(j- 1). Note 
that the bubble is assumed to be small compared to R (oth- 
erwise there can be a significant change in the density of the 
gas as it falls in around the bubble). This would rarely be sig- 
nificant in a cluster, but when it is, then some of the potential 
energy goes into readjustment of the atmosphere as the bubble 
moves. 

For a relativistic gas, 7 = 4/3, so that the enthalpy is 4pV. 
The region where this is dissipated by an adiabatic bubble 
is determined by the pressure distribution of the atmosphere. 
For clusters such as Hydra A and Perseus, roughly half of 
this energy would be dissipated inside the cooling radius. It 
is likely that the bubbles are not entirely adiabatic. On the 
basis of our numbers, radio losses are generally negligible, 
but pieces may be broken away from bubbles, and the rela- 
tivistic particles may leak. Such effects will generally lead 
to a greater proportion of the bubble energy being deposited 
within the cooling radius. 

It is important to note that our estimate of the mechanical 
luminosity relies critically on the assumption that the bubbles 
are close to local pressure equilibrium. This is at least ap - 
proximately true for the Hydra A cluster dNulsen et al.1 2002). 
However, according to the standard view of radio sources, 
bubbles may hav e been significantl y overpressured while be- 
ing formed (e.g. Heinz et al. 1998). In that case, the expand- 
ing bubble drives a shock, and the energy deposited by the 
expansion can be substantially larger than pV . 

There may be additional heat input from the AGNs as- 
sociated with radio outbursts. This could take the form of 
spherical shocks (driven by poorly collimated outflows), di- 
rect injection of relativistic particles, inverse Compton heat- 
ing llCiotti & Ostrikeil 120011) . or other processes. Very sub- 
stantial additional heat inputs would drive convection, leading 
to an isentropic core and mixing out of abundance gradients 
(Bruggen 2002), but this is not a very strong constraint. If 
such energy injection is significantly more than the bubble en- 
ergy input, then it is inappropriate to associate it directly with 
the bubbles, but the mean heating power may be correlated 
with bubble mechanical power. 

Finally, it should be noted that even for adiabatic bubbles, 
the free energy of a bubble decreases with time, and bubbles 
may even break up quickly, so that they disappear as X-ray 
cavities. This means that the instantaneous estimate of bubble 
mechanical power that we have used varies with time, and 
may vary dramatically. A much better controlled sample is 
needed to investigate such issues. 



5.3. Can Cavity Production Quench Cooling Flows? 

We now turn to the question of whether radio sources de- 
posit enough energy into the ICM to quench cooling. We use 
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"X ^spec ^ r cool' 1 
Mechanical luminosity vs. total luminosity minus the spectro- 
scopic estimate of the cooling luminosity. Lines denoting Lx~ Z-spec = Anech 
are shown for the assumption of pV, 4pV, and l6pV energy in the bubbles. 
Symbols and error bars as in FigureQ the arrow denotes an upper limit. 



Lx, the total luminosity of the X-ray-emitting gas from within 
the cooling radius, as an estimate of the classical, or mor- 
phological, cooling luminosity in the absence of heating and 
i-spec, the spectral estimate of the cooling luminosity within 
the cooling radius, as the luminosity of the gas cooling to low 
temperatures. The cooling luminosity, Lx~L spec , must be off- 
set by heating in order to prevent the gas from cooling to low 
temperatures. We note that this quantity ignores non-X-ray 
cooling, such as ultraviolet and optical emission, predicted 
to result from cooling by thermal conduction inside magnetic 
flux loops (Soker 2004), or along reconnected magnetic field 
lines between cold clouds and the ICM iSokeret al] l2004). 
Any such emission would lower the cooling luminosity which 
must be balanced by heating. Figure[2]shows the mechanical 
luminosity plotted against Lx - L spec for our sample. The di- 
agonal lines represent equality between cooling and heating, 
assuming energy inputs of pV, 4-pV, and !6pV per cavity. The 
data are derived from Tables |2] and [3] For RBS 797, an up- 
per limit is shown. The cooling luminosity for RBS 797 is 
poorly constrained by the spectrum, which consists of only 
~ 9000 counts after cleanin. RBS 797, while very luminous, 
is the most distant cluster in our sample and has the shortest 
exposure time (see Section[2j. 

Figure[2]shows several objects, such as Hydra A, Cygnus A, 
and M84, whose cavities can contain enough energy to bal- 
ance radiative losses, at least temporarily, with nearly 1 pV of 
heat input per cavity. The remaining objects, which require 
between a few and ~ 20pV per cavity to balance cooling, 
would do so with varying degrees of difficulty. As discussed 
above, ~ 2pV would be deposited within the cooling radius by 
an adiabatic bubble containing relativistic plasma. Up to 4pV 
is available if the cavities are relativistic and nonadiabatic, and 
there may be further energy input if they are overpressured or 
produce a shock when they are formed. Therefore, the objects 
that require ~ 4pV or less may reasonably be supplied with 
enough energy in the cavities to balance cooling, depending 
on the detailed dynamics (the heat also needs to be distributed 
inside the cooling radius to match the distribution of radiative 



losses). Provided that the true radio cycling timescale ranges 
between f Cs and f r , the cavities in one quarter to one-half of the 
objects in our sample contain enough energy to offset radia- 
tion losses. This would be true of the cavities in the remaining 
obje cts o nly if they are significantly nonadiabatic, as outlined 
in § 15.31 Bear in mind that our conclusions depend on the 
adopted cooling radius (see §3.2), measurement uncertainties 
in the cavity sizes, and the cavity production timescale. Nev- 
ertheless, we can safely conclude that cooling can plausibly 
be balanced by bubble heating in some, but not all, systems. 

It is unnecessary to balance the entire luminosity, Lx~L spec , 
by bubble heating alone if there are other forms of heating 
present. A possible source of heat ing is thermal conduction, 
which, as demonstrated by Voigt & Fabian (2004), could sup- 
ply a significant amount of heat. Using a sample similar 
to our own, Voigt & Fabian found that thermal conduction 
can reduce the cooling luminosity by factors of ^2-3 in 
some objects. Although they are difficult to find in X-ray 
images, shocks associated with the expanding cavities can 
deposit additional energy into the ICM. Deep Chandra im- 
ages of a growing n umber of objects, includi ng Cyg nus A 
( [Wilson et all 120 03). NGC 4636 (Uones etalJl2002l . M87 
(Forman et al. 2004), and Perseus (Fabian et al. 2003a), show 
surface brightness discontinuities that may be associated with 
weak shocks. In Cygnus A and M87, the shocks imply that 
the radio source may provide several times the upper limit of 
the luminosity seen in the bubbles, under the assumption of 
4pV o f energy per bubble (IWilson et alJl2003l iForman et 
120041) . It may therefore be a combination of heating mecha- 
nisms th at leads to quenched cooling, as suggested by several 
authors dBrighenti & Matthewslp002| l2003t iKim & Naravanl 
l2003tlRuszkowski & Begelmanl2002l) . 

It is important to note that that our sample is biased toward 
systems with visible evidence for X-ray cavities, and does not 
represent clusters as a whole. Many clusters, including some 
with l arge cooling flows, do not contain cavitie s (e.g. Abell 
1068, IWise et aLllooll iMcNamara eTai] 120041) . These ob- 
jects may have very different reheating histories than the ob- 
jects discussed here. In this sense, the objects presented here 
represent the best-case examples for reheating the ICM by en- 
ergetic bubbles. Our analysis does not imply that all cooling 
flows can be quenched in this fashion. 

5.4. Trends between X-ray and Mechanical Luminosities 

Figure 13 shows a trend between the X-ray luminosity and 
bubble mechanical luminosity, with the sense that systems 
with larger X-ray luminosities also have larger mechanical 
luminosities. This trend extends over a dynamic range of 
~ 1000 in both X-ray and mechanical luminosity. Just such 
a trend would be expected were the cooling and heating 
of the ICM coupled in some fashion. Several studies (e.g. 
[ Rosner & Tucker! | 1989[ iBinnev & Taborl 119951 iDavid et al] 
l2001tlOuiris et all2001t lChurazov et al. 2002) have proposed 
that cooling is balanced by heating in a self-regulated feed- 
back loop. The feedback loop is driven by episodic radio 
activity fueled by cooling and accretion onto a central black 
hole. The accretion energy is then returned to the ICM 
through an AGN outburst, including the action of the radio 
cavities, which temporarily arrests cooling. At later times, 
the center of the system settles down and the cooling flow is 
reest ablished. During the cooling cycle, mol ecular ga s ( Edge 
2001 ) accumulates and star formation ensues (Johnston e et alJ 
Il987t IMcNama ra & 0'Connelilll989ft . albeit at subst antially 
lower levels than expected in steady-cooling models (Fabian 
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1994). Even if the radio bubbles are not the main source of 
heat, Figure |2] suggests that AGN feedback is intimately in- 
volved in the process that prevents a large cooling flow from 
forming. 

The apparent correlation in Figure|2]should be treated with 
caution. As noted earlier, our sample was selected from clus- 
ters in the Chandra archive with fairly obvious cavities in their 
cores, and neglects those without obvious cavities. Other clus- 
ters are known to have substantial cooling luminosities com- 
mensurate with the observed levels of cold gas and star for- 
mation, yet contain no cavities and have low radio power. A 
prime example is A1068 (Wise et al. 2004; McNamara et al.l 
2004). Similar objects would appear in the lower right of 
this diagram, tending to weaken the correlation. On the other 
hand, it is unlikely that we would have missed objects with 
powerful cavities, which would lie in the upper part of this di- 
agram. Therefore, the distribution of points may represent an 
upper envelope in mechanical luminosity as a function of X- 
ray luminosity. Such a distribution would be consistent with 
the feedback hypothesis, if objects like A1068 are in an ex- 
tended cooling phase in which the central galaxies have ex- 
perienced substantial levels of accretion in the past 100 Myr 
or so, when the radio source has not had a chance to create 
cavities capable of reducing or quenching cooling. 

We have investigated the degree to which other systematic 
effects may lead to an unphysical l uminosity-luminosity cor- 
relation. For example, lElvis et ail J1978I) pointed out that a 
sample of objects with a small range of fluxes and a large 
range of distances will show a correlation in a luminosity- 
versus-luminosity plot, even if there is no intrinsic correlation 
in the sample. Our sample, however, has a large range of ra- 
dio fluxes (from ~ 10 mJy to ~ 10 6 mJy). The cavities in 
our sample also cover a large range of projected angular sizes 
(from ~ 3" to ~ 35"). We believe, then, that these potential 
effects are unlikely to account entirely for the trends seen in 
Figures ^and 13 

Selection bias may also contribute to the correlations. 
Small cavities are easily overlooked in distant objects, since 
cavities of a given linear size become more difficult to de- 
tect as their angular sizes decrease with increasing distance. 
Conversely, in very nearby objects, such as M87, we may 
miss larger bubbles that lie outside the detector. Furthermore, 
other considera tions, such as the bubble position, affect the 
detectability (Enfil in & Heinzl 120021) . The consequences of 
these and other effects on our selection function will be ad- 
dressed in the future using a larger and better-defined sam- 
ple of clusters, including a more sophisticated approach to 
placing limits on cavities that may exist in clusters but were 



missed by the observations. 



5.5. Summary 

We have presented an analysis of 18 systems taken from the 
Chandra archive having clear evidence for cavities in their X- 
ray emission. We find that the energy associated with the cav- 
ities is sufficient to substantially reduce or quench cooling in 
nearly half of the objects in our sample. However, this mech- 
anism alone probably does not provide a general solution to 
the cooling problem, unless X-ray cavities probe only a small 
fraction of the total kinetic luminosity of radio sources. In ad- 
dition, we have discovered a trend between the cooling X-ray 
luminosity and the mechanical energy of the cavities, with the 
sense that more luminous systems produce larger and more 
energetic cavities. The trend, or envelope, may have been 
established by a self-regulated cooling and feedback mecha- 
nism acting in many systems. The existence of such a mecha- 
nism in relatively nearby clusters, where the detailed physics 
can be examined, may provide significant insight on the pro- 
cess of galaxy forma tion that prevails at large redshifts (e.g. 
iVoit & Ponmanl 120031) . A similar mechanism may regulate 
the growth of galax y halos during the dissipative stages of 
their development (Dubinski 1994) and may be an agent re- 
sponsible for the detailed correlation between black hole mass 
and velocity dispersion of spheroids ( Fabian et al. 2002bj). We 
have measured for the first time the distribution of the ratio of 
kinetic luminosity to monochromatic radio luminosity for a 
sample of radio sources. The ratio varies widely, with most 
objects ranging between few and a few hundred, assuming 
IpV of energy per cavity. X-ray cavities provide a unique 
probe of the mechanical power of radio jets, independently of 
the radio properties themselves. 

Our future plans include expanding the sample size and ac- 
quiring better and more uniform radio data. In addition, we 
plan to extend our understanding of the detectability function 
of bubbles, using simulations of images with a wide range of 
exposure and signal-to-noise ratios. 
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sion. This work was supported by Long Term Space Astro- 
physics Grant NAG5-11025, Chandra Archive Grant AR2- 
3007X, and a grant from the Department of Energy through 
the Los Alamos National Laboratory. 
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